SIMPLE UNIVERSAL BOUNDS FOR CHEBYSHEV-TYPE 

QUADRATURES 



RON PELED 

Abstract. A Chebyshev-type quadrature for a probability measure a is a dis- 
tribution which is uniform on n points and has the same first k moments as a. 
We give an upper bound for the minimal n required to achieve a given degree 
k, for a supported on an interval. In contrast to previous results of this type, 
our bound uses only simple properties of a and is applicable in wide generality. 
We also obtain a lower bound for the required number of nodes which only 
uses estimates on the moments of a. Examples illustrating the sharpness of our 
bounds are given. As a corollary of our results, we obtain an apparently new 
result on the Gaussian quadrature. 

In addition, we suggest another approach to bounding the minimal number of 
nodes required in a Chebyshev-type quadrature, utilizing a random choice of the 
nodes, and propose the challenge of analyzing its performance. A preliminary 
result in this direction is proved for the uniform measure on the cube. Finally, 
we apply our bounds to the construction of point sets on the sphere and cylinder 
which form local approximate Chebyshev-type quadratures. These results were 
needed recently in the context of understanding how well can a Poisson process 
approximate certain continuous distributions. The paper concludes with a list 
of open questions. 



1. Introduction 

A quadrature formula is a way of approximating a distribution by a set of point 
masses which preserves the integral of all polynomials up to a certain degree. More 
precisely, given an integer k > 1 and a measure o on R with finite first k moments, 
a quadrature formula of (algebraic) degree at least k is a set of nodes {xj}" =1 C R 
and weights {mj}" =1 C R + such that 



/n 
x j da(x) = 



mix] 



for all integer < j < k. The degree is exactly k if equality does not hold when 
j = k + 1. Such formulas have many applications in numerical analysis, classical 
analysis [15J, geometry [13J and other fields. The maximal degree possible for a 
quadrature formula with n nodes is 2n — 1 (unless a itself is atomic with n nodes 
or less). This degree is attained uniquely for a distinguished formula called the 
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Gaussian quadrature formula whose n nodes are placed at the roots of the n'th 
orthogonal polynomial corresponding to a. 

In this paper, we are concerned with a restricted class of quadrature formulas. 
We consider only probability measures a and restrict our formula to have all its 
weights equal (to -). Hence our formula takes the form 



where the nodes (a?i)£=i nee d not be distinct. Such formulas are called Chebyshev- 
type quadrature formulas. The special case when n = k is known as Chebyshev 
quadratures, see the survey [llj. There is also a multidimensional analogue of ([2]), 
called Chebyshev-type cubatures, when a is a measure on R d , (xj)" =1 C M d and we 
require that f p(x)da(x) = ~ Yli=iP( x i) f° r a ^ polynomials p of degree at most k. 
These formulas arise in various applications, such as combinatorics, statistics [23] , 
potential theory and geometry [28J. In addition, they recently proved essential to 
the understanding of fine properties of the gravitational allocation [61 [7] , where it 
was necessary to understand how well, and with what probability, can a Poisson 
process approximate a given continuous distribution. 

The following questions arise naturally: Does a Chebyshev-type quadrature 
always exist for given a and kl How many nodes are required to achieve a given 
degree for such formulas? 

Definition 1.1. For a probability measure a on K and integer k > 1, define n° a (k) 
to be the minimal number of nodes n required in a Chebyshev-type quadrature ([2]) 
of algebraic degree at least k, or oo if no such quadrature exists. Define n a {k) to 
be the minimal integer such that for any n > n a (k) there exists a Chebyshev-type 
quadrature ([2]) of algebraic degree at least k having exactly n nodes, or oo if no 
such integer exists (see Theorem \l.l\ below). 

Of course, we always have n° (k) < n a (k) < oo (see Theorem 1 1 . 1 1 for an example 
where they have different orders of magnitude). 

The existence question for Chebyshev-type quadratures has been researched ex- 
tensively and is well understood (see [2H [22l [16] ) . Results exist for more general 
formulas than (T5]), involving more general spaces than R and more general func- 
tions than xK In the case of ([2]), one has the following necessary and sufficient 
conditions. 

Theorem 1.1. Given an integer k > 1 and a probability measure a on K with 
J \x\ k da(x) < oo. 

(1) If a is purely atomic with m atoms and k > 2m then the only quadrature 
formula ([I]) of degree at least k for a is a itself. Thus, in this case, if a 
has an atom of irrational weight then n° a (k) = n a (k) = oo and if all atoms 
of a have rational weights then n° a {k) < oo and n a (k) = oo. 

(2) If a either has a non-atomic component or it is purely atomic with m atoms 
and k < 2m, then n a (k) < oo. Furthermore, in this case, there exists an 
n G N such that for any n > n there exists a Chebyshev-type quadrature 
formula t^j for a of degree at least k having all distinct nodes. 





for all < j < k, 
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Although not stated explicitly, the theorem follows readily from results of Kui- 
jlaars [IS] combined with classical results in the theory of the moment problem 
|15j . We prove it in Section I2TT1 

Remark 1.1. If the support of a is contained in some interval [a,b] then it is 
sometimes desirable to have a Chebyshev-type quadrature with all nodes distinct 
and in the open interval (a, b) (see, e.g., [H]j. It is also possible to write necessary 
and sufficient conditions for this case, see Remark \2.1\ for details. 

Theorem 11.11 does not address the quantitative question of the dependence of 
n®(k) and n a (k) on a and k, but part 1 of it already shows that unlike the case 
of the ordinary quadrature ([1]), there is no universal upper bound on n° a {k) given 
only in terms of k. Bernstein was the first to discover the surprising fact that even 
for very simple a, n° a may grow super- linearly. In two papers from 1937 jH [2], he 
proves the following result. 

Theorem 1.2. (Bernstein) Let a be the uniform distribution on [0, 1]. Then there 
exist C, c > such that for all k > 1, 

ck 2 < n° a {k) < Ck 2 . 

Aside from Bernstein's result, the asymptotic behavior of n° a (or n a ) has been 
determined in only a few cases; most notably in [17], where it was generalized 
to a subset of the Jacobi weight functions, and in [18], where it was found for 
measures of the form da = w(x)(l — x 2 )~ l / 2 l x( zy_i^dx for w positive and analytic 
on [—1,1]. We mention briefly that some results exist also for Chebyshev-type 
cubatures. There, research has mostly concentrated on the case where a is the 
area or volume measure of a certain set. See [H] (and [19] for related ideas) for 
results on simple two and three dimensional shapes, and [H [5] for recent progress 
on spherical designs, the case when a is the uniform measure on a sphere, a long 
standing open problem. 

There also exist results in the literature: [21] (inspired by [22]) and [26J, giving 
upper bounds on n a (k) for general measures a in some class. However, these 
results require specific bounds on a which seem difficult to obtain for general 
measures. For example, the result of [21] requires, as one of its ingredients, a 
lower bound on the smallest eigenvalue of the matrix A = (aij)ij=i) where : = 
f(x z — rrii){x^ — rrij)da(x) and rrii := J x l da(x). Moreover, the results require a 
to have a certain regularity: to be non-atomic with full support on some interval 
[21] . or to have a density satisfying certain upper and lower bounds [26J. 

This paper has several goals: First, to give an upper bound on n a (k) which is 
given in terms of simple properties of a (Theorems 11.31 H-41 and II. 5p . requiring 
only an estimate on a's inverse modulus of continuity. Moreover, while the bound 
is particularly simple for absolutely continuous measures with bounded densities, 
it extends also to singular measures and even to purely atomic measures, provided 
some control over the size of the atoms is known. We also give a lower bound 
on n° a (k) which only requires estimates on the k — l'st and fc'th moments of a 
(Theorem 11.81) . Corresponding examples illustrate the sharpness of our bounds 
( Theorems 1 1 . 91 and 1 1 . 1 ft . In particular, we find that for measures a supported on 
[0, 1] with essentially bounded density, n a {k) may rise at most exponentially with 
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k, and this rate of growth is attained for some a. As one corollary of our theorems, 
we obtain an apparently new result on the Gaussian quadrature (Corollary II. 7p . 

Second, to introduce the concept of random Chebyshev-type quadratures (and its 
higher-dimensional analogues), where nodes are chosen by independent samples 
from a (Section 11.21) . We explain how this concept provides another way to upper 
bound n° a (k) and n a (k) and propose the challenge of analyzing its performance. 
A preliminary result in this direction is proven, for the case that a is the uni- 
form measure on a cube (Theorem 11.121) . Our analysis proceeds via a local limit 
theorem. 

Third, to describe applications of our theorems to the construction of point sets 
on spheres and cylinders which are local approximate Chebyshev-type cubatures, 
meaning that one may partition the sphere or cylinder to small diameter sets on 
which the point sets are approximate Chebyshev-type cubatures (Theorems 11.131 
and ll.141) . These constructions and the result for the uniform measure on the cube 
mentioned above, were needed recently in the study of the gravitational allocation 
[7] in the context of understanding how well a Poisson process can approximate 
certain continuous distributions. 

These goals are developed in the next three subsections, without proofs. Sec- 
tion [2] contains proofs and supplements. Section [3] presents open questions. 

1.1. Simple bounds for the number of nodes. In this section, we present an 
upper bound on n a (k) which is calculated in terms of simple properties of a. We 
also give a lower bound on n° a {k) which only requires estimates on the k — l'st 
and fc'th moment of a, and give examples illustrating the sharpness of our bounds. 
The information about a we shall need for our upper bound is contained in the 
following function, 



defined for < 5 < 1. R a is the inverse modulus of continuity of <r; R a (S) measures 
the minimal length an interval needs to have in order to have probability at least 
5. 

Theorem 1.3. Let a be a probability measure with cx([0,l]) = 1. Fix an integer 
k > 2 and let 



(3) 



R a (S) : = min (\x - y\ \ x, y e R, a([x, y]) > 8) 




P (_P_ 



) 



k-l 



r 



6(k + 3) U2e. 
and each p G M fc 



Then for each integer n > r 



-i 



satisfying 





there exist (not necessarily distinct) (xj)" =1 C [0, 1] satisfying 




x i = Pj f or a M integer 1 < j < k. 



i=i 
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The theorem states that if the number of nodes n is large enough with respect to 
k and the quantity R a (jfzg) , then there exists a Chebyshev-type quadrature ([2]) 
having the same first k moments as a. Moreover, for each small perturbation of 
the moments of a, there exists a Chebyshev-type quadrature with these perturbed 
moments. The theorem gives explicit bounds on n and on the size of the allowed 
perturbation. Note that to have a non-trivial bound, we must have R (7 (^^) > 0, 
which is equivalent to saying that a has no atom with mass at least For 
generalizations of the theorem to the case of distinct nodes in (0, 1) and to the 
case of functions other than x\ see Remark 12.21 

Of course, the most important case of the theorem is when the moments of a 
are unperturbed. In addition, in many applications, one is interested in absolutely 
continuous distributions with bounded densities. If the density bound is M, we 
have Rcr(5) > M _1 <5 for all 5. Similarly, if one considers singular a, a typical 
scenario is when R a {$) > cS 13 for some j3 > I. For these cases we have the 
following useful corollary which follows immediately from Theorem 11.31 

Theorem 1.4. Let a be a probability measure with cr([0, 1]) = 1. 

(1) Suppose that a is absolutely continuous with a density which is essentially 
bounded by M. Then for each integer k > 2 we have 

n a {k) < [75e 4 A;M(12eM) (fe ~ 1) l. 

(2) Suppose that R a (5) > c5^ for some c > 0, (3 > 1 and all < 5 < 1. Then 
for each integer k > 2 we have 

n a (k) < [a fc (A; + 3) (/3 - 1)fc+1 l, 

where a > depends only on c and f3. 
Furthermore, in both cases we have that all quadrature nodes lie in [0, 1]. 

Hence, for measures with bounded densities, one needs at most an exponential 
number of nodes in a Chebyshev-type quadrature. A more singular measure may 
require even more nodes. Similar conclusions can be drawn for other measures 
according to which lower bound one has for R a . 

The previous theorems provide quantitative bounds for n a (k) in the cases when 
a does not have large atoms. Can we provide similar bounds when a is a mixture 
of a large atom and a non-atomic component? or when a has infinitely many 
atoms? The following theorem does so. Define, for a probability measure a and 
< e < 1, 

°~ (<r({x}) - e)6 x , 

{x I a({x})>e} 

In words, cr* is a with all its atoms truncated to mass e and a' £ is its normalized 
version. 



< ■= 
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Theorem 1.5. Let a be a probability measure with cx([0, 1]) = 1. Fix an integer 
k > 2 and suppose that there exists < e < 1 such that 

/_\ e 2 

( 5 ) ,mn m < 



Fir siic/i an £ and Zet 



a*([0,l]) 2A; + 7' 



p := {k — l)R a , 



r :- 



P (J-) 
6{k + 3) \12eJ 

1 2k±§\ 



2k + 7 

k-l 



Then for any integer n > max ^ ro .tQ ij) ? ^r^j there exist (not necessarily distinct) 
(xi)i=i Q [0, 1] satisfying 



n 

n ^ 



x{ = J x-'da(x) for all integer 1 < j < k. 



Remark 1.2. (1) It is not difficult to see that condition (jSJ is satisfied for 
any small enough e > if a has a non-atomic component or at least k + 4 
atoms. 

(2) Note that the largest atom in a' e has at most ^ mass so that condition 
(JSJ) ensures that p > 0. T/ie reason that the o/ Theorem \1.3\ is replaced 
by one? /or t/ie extra factor m t/ie bound on n is that we may not 
be able to exactly truncate the atoms of a to probability e using atoms of 
size -. 

n 

(3) Similarly to Theorem ! 1.31 we can quantify a statement saying that for any 
moment vector p which is close enough to the moments of a, we can find 
{ x i\7=i with these moments. 

(4) The proof is based on writing a = qo~\ + (1 — g)o~2 f or probability measures 
o~i,o~2, where o\ approximates a' £ and o~2 is the "leftovers" of the large 
atoms of a. The approximation is chosen so that o~2 already has atoms 
with rational probability, then Theorem M .3\ is used to get a Chebyshev-type 
quadrature for a±. We note that this approach might yield better bounds 
than those of Theorem M . 3\ even for a which do not have large atoms. For 
example, if R a {-^) is very small, one may try to decompose a = gcri + (l — 
g)<72 so that R ai {-^) > R (T (jt^) and q is rational with small denominator. 
Then approximate 02 in a simple manner, say as in Lemma 2.4 below, 
and finally approximate o~\ using Theorem M. 3\ and use the freedom in the 
moments afforded by fll]) to compensate for the errors in the moments of 
the approximation to o~2- 

Lower bounds. In this section we complement the above upper bounds for n a {k) 
by presenting lower bounds for n° a (k) and examples illustrating the sharpness of 
our bounds. As a by-product of our results, we note an apparently new inequality 
for the Gaussian quadrature. 

We start by describing a lower bound for n®(k), for general probability measures 
a, which Bernstein used in deriving Theorem 11.21 To state it, we first recall that 
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for k = 2m— 1 for m G N and a probability measure aonR with J \x\ h da(x) < oo, 
unless a is purely atomic with less than m atoms, we have the Gaussian quadrature 
formula with nodes £^ < £i> < • • • < $}m^ and weights (A- m ' ) )™ =1 satisfying 

m „ 

(6) \ M (d m) ) j = / x 3 da(x) V < j < k. 

i=l ^ 

Theorem 1.6. (Bernstein pQj For a probability measure a and k = 2m — 1 as 

above, we have n° a (k) > . A (m) - ■ 

Bernstein proved this theorem in the special case of the uniform distribution 
on an interval, but as some authors note P2J EH EH], the bound extends to all 
measures. We note an immediate corollary of Theorem 11.41 and Theorem 11.61 to 
an estimate on Gaussian quadratures. 

Corollary 1.7. For any probability measure a with er([0, 1]) = 1 and density 
essentially bounded by M, we have for any m G N that 

Ai m) > 1 



[75e 4 (2m - 1)M (12eM) (2m " 2) ] ' 

This estimate appears to be new and we do not know if it is simple to prove 
directly. Similar corollaries can be phrased for general measures using R a and 
Theorem 11.31 

Theorem 11.61 can be quite accurate (as Theorem 11.21 illustrates), however, one 
drawback of it is that it may be difficult to apply in specific cases since it requires 
knowledge of the Gaussian quadrature associated to the given measure. We now 
propose a second lower bound, whose proof makes use of similar ideas to that of 
Theorem 11.61 which has the advantage that in order to apply it, the only required 
information about the measure are bounds on its (k — l)'st and fc'th moments. 

Theorem 1.8. Let a be a probability measure on R with o"({0}) < 1. Then for 
every odd integer k > 3 for which J \x\ k da(x) < oo, we have 



k-l 



(7) n° a (k) > 



(/ x k da(x)) 
(J x k ~ l da(x)y 



We remark that the lower bound given by the above theorem changes, in general, 
when replacing a by a translate of it. Hence, one may wish to optimize the amount 
by which to translate a before applying the bound. To keep the theorem as simple 
as possible, we avoid making this optimization here. 

Theorem 11.81 implies that, for example, to obtain a lower bound on n®(k) for 
probability measures supported on [0, 1] (which are not 5 ), it is sufficient to have 
a lower bound on the absolute value of the fc'th moment of a and an upper bound 
on the (k — l)'st moment of a. The following corollary is proved via this technique. 
It illustrates that the upper bound given by Theorem 11.41 and the lower bound 
given by Theorem 11.81 may be close in specific examples. 
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Corollary 1.9. There exists C > such that for every odd integer k > C, there 
exists a probability measure o~%. on [0, 1], absolutely continuous with density essen- 
tially bounded above by Ck, satisfying 



In fact, the measure constructed in this corollary is simply the exponential 
distribution, properly truncated and rescaled. Note also that, since n° a {k) is non- 
decreasing in k for any measure a, the corollary implies a similar bound for even 
integers k. 

Let us compare the lower bound of Corollary 1 1 . 91 with the upper bound on n a (k) 
given by Theorem II .41 Since the density of is bounded by Ck for some C > 0, 
Theorem 11.41 gives 



for some C > 0, which differs from the bound of Corollary 11.91 by a log k factor 
in the exponent. Seeking to have an example on which the upper bound of Theo- 
rem [T3] is sharp, up to the constants involved, we introduce the following second 
example. We set d n (a), for n e N and a probability measure a on K with all mo- 
ments finite, to be the maximal possible degree of accuracy for a Chebyshev-type 
quadrature for a having exactly n (not necessarily distinct) nodes (or oo if any 
degree of accuracy can be attained). 

Theorem 1.10. Let o" be the probability measure having density 



Thus the theorem shows the surprising fact that for the uniform distribution 
on two disjoint intervals (Tq, d n (o~o) has completely different orders of magnitude 
for odd and even n. It also shows that the upper bound of Theorem 11.41 can be 
attained, up to the constants involved. Our proof of this theorem uses general 
theorems of Peherstorfer [20] and Forster & Ostermeyer [9] which, when taken 
together, show that d n (a) may rise at most logarithmically in n for odd n, whenever 
a is a symmetric measure having outside its support. We remark also that the 
phenomena that d n (a) may have very different orders of magnitude for odd and 
even n, was first discovered, in a particular case, by Forster [8]. 




nUk)<n ah (k)<{Ck) h 
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1.2. Random Chebyshev-type cubatures. We give the name random 
Chebyshev-type cubature to the situation in which we would like to approximate 
the moments of a measure a by the moments of a uniform distribution on n points 
(as in ordinary Chebyshev-type cubatures), but do not choose the position of the 
points, instead, the points are chosen randomly according to independent samples 
from a. In such a situation, it is natural to ask how small is the probability that 
the moments of the random measure approximate the moments of a very well. In 
general, this is a question about a small ball probability. 

As we shall see, this notion gives another way to prove existence of Chebyshev- 
type quadratures and cubatures and we believe that it deserves better study. In 
addition, in the analysis of [7J, it arose naturally in the context of understanding 
how well a Poisson process approximates Lebesgue measure. 

To formalize the above, fix d > 1 and define for k > 1, PolyDim(/c, d) := 
{ k+ d d ) - 1. Then define the moment map P fc d : R d -> RPdyDMM by 

(12) pf(x) : = (x a ) a 

where a runs over all multi-indices with < \a\ < k, where we mean that a G 
(N U {0}) d , x a := nil x? and \a\ := £ii a t . 

Given a probability measure a with support in IR d , let M k (a) denote its vector 
of multi-moments of degree at most k. That is, 

M k (a) := ( / x a do(x) 

' o<H<fc 

For n > 1, consider the following random measure 



1 - 
a n := - y^S. 



n . 

i=i 



where the {xj}" =1 are chosen independently from the distribution a. Note that 
M k (a n ) = ~Y™=i P k( x i)- lt follows that EM k (a n ) = M k {a). Still, the moments 
of o~ n typically do not approximate well the moments of a. Indeed, by the central 
limit theorem, the difference \M k (a n ) a — M k (a) a \ scales like for any fixed a (if 
a has moments of any order, say). We are interested in the probability that this 
difference is much smaller. More precisely, let 

Pn,kM ■= P ( ll M fcK) " Mjfe^Hoo < 4= 

V \/n 



This is the small ball probability for the random vector M k (a n ). We would like 
to understand how it scales for a fixed n as e tends to 0. The following lemma 
connects this probability to the existence of Chebyshev-type cubatures. 

Lemma 1.11. If for some n, k > 1 and every e > we have Pn,k,e{°~) > then 
there exists a Chebyshev-type cubature for a of degree at least k having exactly n 
(not necessarily distinct) nodes. 

Thus understanding p n ^, e provides a different way to show existence of 
Chebyshev-type cubatures (and to prove lower bounds for n°). We propose the 
challenge of bounding, in specific examples, the minimal n for which the condition 
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of the lemma is satisfied and seeing if this approach may improve known bounds. 
The most interesting case in this respect is that of spherical designs, when a is the 
uniform measure on the sphere § d_1 , but one may start by checking what bound 
is achieved for the interval and comparing it with Theorem 11.21 In this paper, 
we content ourselves with a small step in this direction (which, however, already 
takes some work to prove) by showing that the condition of the lemma is satisfied, 
for large enough n, when o is the uniform measure on the cube [—1, l] d . This is 
achieved by showing that, for large enough n, the random vector Mk(<J n ) has a 
positive density at M k (a). Unfortunately, our result does not provide quantitative 
bounds for n. 

Theorem 1.12. Fix k > 1 and let {Xi)^ be an IID sequence of RV's uniform 
on [-1, Let Mi := P fc d (X;) and S n := Y%=i( M i ~ EMi). Then there exists 
Nq = No(k, d) > 0, a = a(k, d) > and t = t(k, d) > such that for all n > No, 
S n is absolutely continuous with respect to Lebesgue measure in R Pol y Dim ( fc ' d ) and 
its density f n (x) satisfies f n (x) > a for \x\ < t. 

This theorem can be viewed as a local limit theorem (and its proof follows this 
approach). The moment vector ^JnS n is a sum of IID contributions and we show 
that starting at some large n, it has density which (suitably scaled) converges 
uniformly to the density of a centered Gaussian vector. Our main tool is Fourier 
analytic estimates. 

1.3. Application to construction of local cubatures. In this section, we ap- 
ply the results of Section 11.11 to give a construction of discrete point sets on the 
sphere and the cylinder which are approximate Chebyshev-type cubatures. In the 
case of the sphere, we approximate its surface measure. In the case of the cylinder, 
the approximation is to a measure with density (with respect to surface area) con- 
stant on every spherical section and growing linearly along the axis of the cylinder. 
In both cases, our approximations are stronger than ordinary Chebyshev-type cu- 
batures in that they are "local", i.e., there is a partition of the set in question 
(the sphere or the mid-part of the cylinder) to subsets of small diameter such that 
our point set restricted to each of these subsets is an approximate Chebyshev-type 
cubature. 

The application to the cylinder, which builds on the application to the sphere, 
was central in the recent study [7J on gravitational allocation where it was used 
to construct "wormholes" ; long tentacles in space surrounded by rings of stars in 
which the gravitational force is atypically strong in the tentacle's direction. 

Our construction of the Chebyshev-type cubature for the sphere is very similar 
to a construction of Wagner [27] which he used in his work on a problem in 
potential theory (in a somewhat similar manner as the application in [7J). Despite 
this similarity, we chose to give a full proof of it here since some parts in Wagner's 
construction (such as the exact partition of the sphere) are only sketched and since 
our construction gives explicit bounds on the number of nodes in the cubature (at 
the expense of getting only an approximate cubature formula), whereas his only 
shows existence. 
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To state our theorems, define a d to be the d- dimensional area measure on sets 
in R d , and, abusing notation slightly, also as the d- dimensional Hausdorff measure 
on sets in R d for d' > d. For a set E C M. d , define Diam(.E') := max x y ^E \ x ~ y\i 
the diameter of E, where | • | is Euclidean distance. We use the same multi-index 
notation as defined after (fl2|) . Define for d, k > 1 and 5 > 0, 

/ fee \ m+1 A 
(13) m (d,k,5) := Smallest integer m > 1 satisfying I J < 



m + 1 J ~ 2d2 k ' 
Embedding E> d into M. d+1 as the unit sphere, we prove: 

Theorem 1.13. For each d > 1, k > 1 and < r, <5 < 1 £/iere exist C(d), c(d) > 
(depending only on d), an integer K = K(r,d) > and a partition of E> d (up 
to surface measure 0) into measurable subsets E\, . . . , E^ satisfying the following 
properties: 

(1) Diam(Ei) < C(d)r, a d {Ei) > c(d)r d for all i, and K < C(d)r~ d . 

(2) For N = n d , where n can be any integer satisfying n > C (d) mo<yd ' k,5 \ and 

■;j)j = l 

N 



for each 1 < i < K, there exist (zij)f =l C Ei such that 



(14) 



fE^')-^)/ e , 



g(z)da d (z) 



< S 



for all g : M d+1 — > M of the form g(z) = (z — w) a for w G E> d and a 
multi-index a with \a\ < k. 



To state our theorem for the cylinder, we make a few more definitions. Given 
L, W > and a dimension d > 1, let 

Pl,w ■= {xER d | <L,x 2 2 + --- + x 2 d = W 2 }. 

so that Pl,w is the curved part of the boundary of a length L cylinder of radius 
W. Let 

vl,w be the measure supported on Pl,w and absolutely continuous with 
respect to cr^-i with density V(xi, . . . ,Xd) = v(xx) = 1 + S2 ^ ± - he., the density 
increases linearly from 1 to 2 as x\ increases from — L to L. Recalling the definition 
of mo(d, k, 5) from (TIB"]) , we prove: 

Theorem 1.14. For each d > 3 there exists C(d) > such that for each k > 
1, L > C(d), W > 0, < t < W and < 5 < W k , we have an integer 
K = K(L, W, t, d) > and measurable subsets Di, . . . , Dk Q P2L,w satisfying the 
following properties: 

(I) V2L.w{Di n Dj) = for each i ^ j and v 2 l,w {Pl,w \ (u^A)) = 0. 
(II) Diam(Di) < C{d)r, v 2 L,w{Di) = r^ 1 for alii, and K < C(d)LW d - 2 T^ d - l \ 
(III) Forn = n d ~ x , whereni can be any integer satisfying n% > C(d) mo ^ d ~ 2 ' k ' 5 ^ 4LWS) \ 
and for each 1 < i < K , there exist (wD it j)j = i ^ A such that 



1 *. n ^ 1 /" 

-yiKwDij) 77TT / K w ) d ^L,w( w ) 

n J^! v 2 L,w(Di) Jd, 



< S 



for all h : M d — > E of the form h{w) = {w — y) a for y G P 2 l,w an d a 
multi-index a with \a\ < k. 
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It is worth noting that in the above theorem, the sets Di, . . . , Dk cover Pl,w 
(up to v- measure 0) and are contained in P2L,w-, but do not necessarily form a 
partition of Pil,w ( u P to //-measure 0). Indeed, this is not possible for generic 
values of L, W and r since property (I) implies the Di are disjoint (up to //-measure 
0) and property (II) implies that each Di has u measure exactly r d ~ x . Thus, for 
the sets to form a partition, we would need that V2L ' W J£^' W ^ be an integer. 

1.4. Acknowledgments. We are grateful to Franz Peherstorfer for explaining to 
us the relevance of [20] and suggesting that it may be used to obtain lower bounds 
for n a (k), as we did in Theorem 11.101 We thank Nir Lev for his essential help in 
referring us to the book of Stein and explaining the relevance of its propositions 
to estimates in Section 12.31 We thank Boris Tsirelson and Mikhail Sodin for 
several useful conversations, in particular concerning approximation of continuous 
measures by discrete ones. Finally, we thank Greg Kuperberg, Yuval Peres, Dan 
Romik and Sasha Sodin for useful discussions and comments on quadratures and 
cubatures. 

2. Proofs and supplements 

2.1. Existence of Chebyshev-type quadratures. In this section we prove 
Theorem 11.11 We start with the following theorem which follows from classical 
results in the theory of the moment problem: 

Theorem 2.1. Given k = 2m — 1 for m G N and a probability measure a on R 
with J \x\ k da(x) < oo. Unless a is purely atomic with less than m atoms, there 
exists a quadrature formula §\§ for a of degree at least k having exactly 2m + 1 
nodes. 

Proof. We assume without loss of generality that a is atomic with m nodes since 
otherwise we can replace a by its Gaussian quadrature ([6]). Then its support is 
contained in an interval [a, b}. Fix c < a and d > b. Note that for any polynomial 
P ^ of degree at most k which is non-negative on [c, d] we have f Pdo~ > 0. In 
other words (see [T31 III §1]), a (or rather its first k moments) is strictly positive 
with respect to [c, d] and k. This implies that there exists a quadrature formula 
o~ m+ \ for a having degree at least k, exactly m + 1 nodes, all in (c, d) and all 
different from those of a (this is any of the lower representations of index n + 3, 
see [T5| III §7.1]. All nodes are in the interior of [c,d] since k is odd). Then the 
measure |(cr + cr m+1 ) satisfies the requirements of the theorem. □ 

We also need two theorems of Kuijlaars: 

Theorem 2.2. (|16j. Theorem 3.2) Given a probability measure a, suppose we 
have a Chebyshev-type quadrature formula ()2]) of degree at least k with n nodes 
in (—1, 1), of which n > k are distinct. Then there exists a Chebyshev-type 
quadrature formula with n distinct nodes in (—1, 1) of degree at least k. 

Theorem 2.3. (|16j. Theorem 4-2) Given a probability measure a and a quadra- 
ture formula ([1]) with weights (mj)" =1 , distinct nodes in (—1, 1) and degree at least 
n — 1, there exists a relatively open subset U of the collection {(pi, . . . ,p n ) I Pi > 
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for 1 < i < n and ^2™ =1 Pi = 1} with (mj)™ =1 G U and such that for every 
(Pi)i=i U there exist nodes (5j)" =1 C (—1, 1) satisfying 



The second theorem was proven in [16J for absolutely continuous a with bounded 
support but the (short) proof is valid for any a. 

Proof of Theorem li.il If a is purely atomic with j ' < | atoms, it is well known 
that a itself is the only quadrature formula having degree at least k. This can be 
seen by considering the non-negative polynomial P having a double zero at each 
atom of a. Since P has degree 2j < k and since J Pda = we see that the integral 
of P is zero also with respect to a quadrature with degree at least k. Hence, the 
nodes of that quadrature are a subset of the nodes of a, but this implies that they 
are equal since the location of the nodes determines the weights by solving a linear 
system with a Vandermonde coefficient matrix. 

Assume now that a is not purely atomic with j < | atoms. We may assume 
k = 2m — 1 for some m G N since if k is even, the theorem remains true when k 
is replaced by k + 1. We use Theorem 12. II to obtain 02m+i, a quadrature for a of 
degree at least k having exactly 2m + 1 nodes. It follows from Theorem 12.31 that 
for some hq G N and any n > hq there exist Chebyshev-type quadratures with 
n nodes having the same first 2m moments as cr 2rri+ i, so that, in particular, they 
have degree at least k with respect to a. The nodes of these quadratures can be 
made distinct using Theorem 12.21 □ 

Remark 2.1. For a probability measure a with support in [a,b] and k G N, we 
say that a is singular with respect to [a, b] and k if there exists a polynomial P ^ 
of degree at most k which is non-negative on [a,b] and such that J Pda = 0. 
Equivalently, o is singular with respect to [a, b] and k if and only if it is purely 
atomic and its index 1(a) < k where 1(a) := ^2 I(x) over all atoms x of a and 



Kx) 



1 x = a or x = b 

2 x G (a, b) 



If a is singular, it follows from the same proof as above, but using the polynomial 
P exhibiting the singularity, that the only quadrature for a with all nodes in [a, b] 
is a itself. 

If a is not singular with respect to [a, b] and k, then the same proof as above with 
minor modifications shows that for any large enough n, there exist Chebyshev-type 
quadratures for a having degree at least k and n distinct nodes in (a, b) . If k is 
odd, the only modification is that in the proof of Theorem \2.1\ one should obtain 
its representations directly in [a, b] without passing to the larger interval [c, d] (this 
is possible since a is non-singular on [a, b}). Then the quadrature thus obtained 
will have all its nodes in (a, b). For even k, the additional required modification is 
to first replace a by a' , a canonical representation of it with support in [a, b], index 
k + 2 and the same first k moments [Tol III §4], then to apply the above proof for 
a' and k + 1. 



14 



RON PELED 



2.2. Bounds for Chebyshev-type quadratures. In this section we prove The- 
orems 11.31 and 11.5} which give general upper bounds for the minimal number of 
nodes required in a Chebyshev-type quadrature, and Theorem 11.8} Corollary 11.91 
and Theorem 11.101 which give lower bounds for the required number of nodes and 
examples of cases where many nodes are required. We remark first about possible 
generalizations of our results. 

Remark 2.2. (1) It is sometimes desirable that the nodes of the quadrature 
be distinct and contained in the open interval (0, 1). To obtain bounds for 
such formulas using our results, start by picking a small e > 0, linearly 
map a to have support in [e, I — e], apply Theorem \1.3\ to the new measure 
and use the freedom afforded by to make the moments of the resulting 
quadrature equal those of a. To make the nodes distinct, use Theorem \2.2\ 
(the proof of Theorem M.'d gives at least k distinct nodes). 
(2) It may also be desirable to have a result similar to Theorem \1.3\ for func- 
tions other than x J . The main ingredient required to adapt our proof to 
such a setting is to have a "quantitative inverse mapping theorem", as in 
Proposition \2. 51 for the new collection of functions. 

We start our proofs by recalling the definition ([3]) of R a for a probability measure 
a, and noting the following simple properties: 

(1) R a is monotonically increasing. 

(2) R a (5) = if and only if a has an atom of mass at least S. 

(3) If a is supported on [0, 1], then R a (~) < — for each m G N. 

(4) If a is absolutely continuous with a density that is essentially bounded by 
M, then R a {8) > £ for all 5. 

2.2.1. Proof of Theorem \1.3[ We start with a lemma providing a simple approxi- 
mate Chebyshev-type quadrature. 

Lemma 2.4. Let a be a probability measure with er([0, 1]) = 1. For n G N, let 

\i :— i Y^7=i^Vi' w here the (y 4 ) are chosen according to the rule: 

(15) yi :=mm(y(= [0,1] | <r([0,y]) > l<i<n. 



Then for all j G N, we have 

x 1 do — I x^djjL 



< 



1 

n 

Proof. Define yo ■— and for each < i < n define the "leftover mass at y" by 

a t := a([0,yi\) - -. 

n 

Note that < a>i < by definition of the y^. Also, define measures (c i )^ =1 

by 

<7i(-) := a(- n {yi-iM + 0^-1^(0 - di5 yi {-) 
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and note that these measures are non-negative with total mass exactly ^ and that 
we have a = J2i=i a i- Now, fix j G N and estimate 

/ x ] do - j x°dp = f{x J - yDdai < - (yf - yj^) < - 
J J i= i J n i= i n 



as required. 



□ 



Remark 2.3. It is worth noting that - for some c > is the best approximation 
possible in this level of generality if one uses the above method of dividing a into 
ai 's with mass - and approximating each one with one point. This can be seen 
by considering the example of a = ^(Sq + S\) when n is odd and the example of 



a 



|(5o + Si + Si) when n is even. 



3^ - 2 

Our aim is to perturb the above simple approximation into a Chebyshev-type 
quadrature for a. To this end, we define the moment map T fc : R k -> R k by 

k 



i=l 



and rely on the following quantitative "inverse mapping theorem" : 
Proposition 2.5. Fix p > 0, integer k > 2 and let z6R fc satisfy 



(16) 



9 



3(k 



< Zi < 1 - 



9 



3(k-l) 



and 



9 



k-1 



for all 1 < i, j < k with j ^ i. Then for any p G M. k satisfying 

k-l 



(17) 

there exists w G 



\P-Tk(z)\ 
i k satisfying \w — z 



< P_(_P_\ 
00 " 3 \12e) 



< E 

00 — 3(fc-l) 



and Tk{w) = p. 



In words, the proposition shows that if the Z{ are well separated, then the image 
through Tfc of a ball around z contains a ball in moment space (where the balls 
are in the l^ metric), and it gives quantitative bounds on the radii of these balls. 
Since the proof of this proposition is somewhat long, we delay it until after we 
explain how Theorem 11.31 follows from the proposition and lemmas. 

Iterating the proposition, we obtain the following corollary. 

Corollary 2.6. Given fi := - J^ILi all < y,i < 1 and k > 2. Suppose 

that there exist < p < 1 and s disjoint subsets (z(r)) s r=1 of the (?/«) 's, each of 
size exactly k, such that 

9 / / 1 P 



3(fc — 1 

\z(r)i - z(r)j\ > 



< z{r) t < 1 
9 



3(fc - 1) 



k 



for all 1 < r < s and 1 < i,j < k with j ^ i. Then for any p G IR fc satisfying 
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there exists p! of the form p! := \ YHi=\ $xi w ^h a ^ < %i < 1, swc/i £/ierf J a^c?// = 
Pj for all j. 

We emphasize that in this corollary and the rest of the proof, by disjoint subsets 
z(r) C (i/i) we mean that we may choose indices . . . , i r k ) s r =i such that z(r)j = y^ 
and each % appears at most once in all these index sets. Note that with this 
convention, if the (yi) contain a certain value multiple times, then it may happen 
that the (z(r)) also contain this value multiple times. 

Proof. The corollary follows by applying Proposition 12.51 to each of the subsets 
z(r), each time changing the moments of the measure in the direction of the 
vector p. Note the additional factor ^ in ( Tl9j) as compared to ffTTj) . This factor 
appears since is an unnormalized sum whereas p! contains the normalization 
factor -. □ 

n 

Finally, it remains to show that if a measure does not have large atoms, then the 
simple approximation of Lemma 12.41 contains many disjoint subsets as in Corol- 
lary E5J 

Lemma 2.7. Let a be a probability measure with cr([0, 1]) = 1. For n G N, let 

(?/j)™ =1 be the simple approximation ([151 . Then for each integer k > 2 such that 
n > k(k+3) there exist [jrgl disjoint subsets (z(r)) of the (yi), each of size exactly 
k, which satisfy ([TBI) with p = (k — l)R a {-^) . 

Proof. Define yo := 0. Note that by definition, we have yj — yi > R,j(—) for 
< i < j < n. Let io := \]^hf\ and define the subsets (z(r)y° =1 by 

^v)j • Vjio+i — 1 

for 1 < j < k. We need to verify that the conditions in ( fl8l) hold with p = 
(k — l)R a (-j^). To check the first condition, note that since R a is non- decreasing 
and ^3 < i < fe^g + 1, we have 

z(r)i > z(l) x = y i0 > R a f - ) > R a ( — j— ) > ^ ^ ' 



nj ~ \k + 3J ~ 3 V fc + 3 

. . . . / (k + l)i - 1\ 
z(r)j < z(z )fc = 2/(fc+i)i -i < 1 - it<r I 1 I < 

fc + l fc\ / 1 \ 1 

< 1 - R a 1 - — — - - < 1 - RA — — < 1 - -ik 



using the assumption that n > kik + 3). The second condition in ( fl8l) follows 
similarly. □ 

Remark 2.4. VFe note that there do exist a with atoms of size for which the 
(Vi) °f p5p do not contain even one subset which satisfies the separation condition 
(JTSJ for a positive p. For example, a = Yli=i Hence the above lemma is 
close to optimal. 

Putting all the above claims together, we may finish the proof of Theorem 11.31 
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Conclusion of the proof of Theorem \1.3l Let p and r be as in the theorem. Fix 
an integer n > r _1 and a vector p G M fc satisfying (J3J). By Lemma 12.44 we have 
(yi)i=i Q [0) 1] sucn that for all j e N we have 



(20) 



n 

t=i 



x^da(x) 



1 

< -. 

n 



Note that using the 3'rd property of appearing in the beginning of the section, 
we have p < 1 and so n > fc(A; + 3). Hence, by Lemma [2.71 there exist s := [^3] 
disjoint subsets (z(r))* =1 of the (?/«), each of size exactly k, which satisfy ( ITS]) for 
the given p. Hence, by Corollary 12.6} for any p' e R fc satisfying 



(21) 



J n *— ' 



i=l 



~ 3n V12e/ 



fc-i 



for all j, there exist (xi)™ =1 C [0,1] such that ^Yli=i x i = P'j f° r a ^ 1 — J ' — k- 
Since p satisfies (TjJ, equations ( 12"U]) and f[2"T]) will imply the theorem if 

k-l 



+ l <£?_(_£_)> 

- 3n \12eJ 



n 



This now follows by the definition of s and the condition n > r . 



□ 



Proof of Proposition Fix an integer k > 2. We first define some notation: for 
w e 



l k , let V(w) be the Vandermonde matrix defined by 

/I •■■ 1 \ 
wi ■■■ w k 



V(w) 



k-l 



W, 



k-l 



and let U (w) be a slightly modified version defined by 

U(w) : 



/ 1 

2wx 



1 \ 

2w k 



k-l 



For a matrix A e M kX k, define ||^4||oo := niaxi<j<jt X^Li l^k?l> the infinity norm of 
the matrix. We continue by citing (a special case of) a theorem of Gautschi about 
norms of inverses of Vandermonde matrices [10] ■ 



Theorem 2.8. ( Gautschi) For w G WL k satisfying w,i > for all i and 7^ Wj for 

all i 7^ j , we have 

k 1 , 



\\V(w)-% 
We immediately deduce 



max 

Ki<k 



n 

i=i 



W ; 
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Corollary 2.9. For w G M. k satisfying Wj > /or a// % and W{ ^ Wj for all i ^ j , 

we have 

h 



l<i<k- L± \Wj — W 

- j= i I * 



In particular, if < Wi < 1 for all i and there exists < a < 1 such that 
\wi — Wj\ > -j^j for all i j then 



II Wll- < I (|) 



fc-i 



Proof. Noting that = -DV^w) where D is a diagonal matrix with 1,2, ... ,k 

on its diagonal, we see that 



II Wi*, < IIvh^iuii^iu = ll^H- 1 ! 



OO ) 



so the first part of the corollary follows from Theorem 12.81 For the second part 
assume for simplicity that k is odd. For the case k — 1 there is nothing to prove 
for k > 3 the assumptions and Stirling's approximation imply 



Ki<k -•- -•- \Wj — W 



:i I 



1 

< 




-i 



n(k - 1) 

Similarly, one can check that the required estimate holds when k > 2 is even. □ 

We continue the proof by defining a vector field G : M, k — >■ M fc and an ODE, 

G(w) := U(w)- X {p - T k (w)), 
w(t) := G(w(t)) and w(0) := z. 

By standard existence theorems for ODEs, there exists a solution to the ODE 
w : [0, r*) defined up to the first time that G(w(t)) is undefined, i.e., the first time 
that Wi(t) = Wj(t) for some i ^ j. r* = oo if such a time does not exist. Let also 
t* be the first time that \w(t) — z\oo = c^riy, or infinity if such a time does not 
exist. It is clear from the separation conditions ffl6|) on the coordinates of z that 
t* < r* with a strict inequality if r* < oo. 

Note that the Jacobian -^T k {w) = U{w). Hence, for each t < r*, 

^-(T k (w(t)) -p) = ^-T k (w(t))w(t) =p- T k (w(t)) 
at aw 

from which it follows that T k (w(t))—p = e - *^^) — p). We deduce that if = oo, 
then since for t < t*, w(t) lies in a compact set, we may extract a subsequence of 
w(t) converging to some w with \w — < ^ k p _^ ■ By continuity of T k , this w 
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satisfies Tk(w) = p as required. Hence, we assume, in order to get a contradiction, 
that < oo. We now calculate 

w(t*) — z = / w(s)ds = / U(w(s))^ 1 (p — Tk(w(s)))ds = 
Jo Jo 

r-t* 

e~ s U{w{s)Y l ds{p-T k {z)). 

Hence, noting that by Corollary 12.91 with o = | we have for s <t* that 

1 /12e\ fc-1 



we obtain (using assumption (JE 



1 /12e\ /" 



1 /12e\ ^ fe-i p p 

fc V p / \12e/ 3 3(fc-l) 
contradicting the definition of £*. Thus the proposition is proven. 

2.2.2. Proof of Theorem \1.5[ Recalling the notation of Theorem 11.51 let us fix 
< e < 1 satisfying 

£ 2 

(22) ^OMl) < ^1 

and 

1 2A; + 6 



(23) n > max , 

l ra*([0,l])' £ 

As noted in Remark 11.21 condition (|22|) implies that the right hand side of (j23|) is 
finite. Let A = {x \ a({x}) > ffrj^} and define a measure 

o"2,n := y~] - [n{cr({x}) - e)J 5 X .. 

In words, o^,™ has an atom for every atom x & A and the mass of this atom is the 
largest multiple of - which is no larger than a({x}) — e. Define also 

Cl,n ■= 0- - a 2 , n - 

Then, by our definitions and f[2"3"j) . we have 

(24) al(B) < a l>n (B) < ^^{B) 

for every Borel set B. We now let q := 0i )n ([0, 1]) so that l — q — o"2,n([0, 1]). Note 
that q > and is a multiple of ^ by the definition of <7 2iri . Letting a' l n := ^p, we 
have 

(25) a = qa[ + a 2n . 
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We claim that there exists a Chebyshev-type quadrature for a' l n of degree at least 
k and having exactly qn (not necessarily distinct) nodes in [0, 1] (qn is an integer!). 
By Theorem II. 3[ we know that such a quadrature exists if 

(26) qn > — 

where 

p' := (k - l)R a[ 

p' ( P 



> 



6(k + 3) \12e 



k + 3 

fc-i 



By (I24I) we have that a' ln ([x,y]) < §j^cr' £ ([x,y\) for any x < y. Hence R a > (5) 
i? CT /(||i|(5) for any < 5 < ~ and in particular p' > p and consequently r' > r 
(p and r were defined in the statement of the theorem). In addition, by ( EH) , we 
have that g > <r|([0, 1]). We conclude that holds by (g3). 

To finish, we have obtained a Chebyshev-type quadrature for o"' l n of degree at 
least k, 

■= — 5^*xj 



for some {xj}^ C [0, 1]. Defining 



A* : = °"2,n + - 

n < J 



n 

i=i 



it is straightforward to check using (125|1 that /x is a Chebyshev-type quadrature of 
degree at least k for a having exactly n (not necessarily distinct) nodes in [0, 1]. 

2.2.3. Lower bounds for the number of nodes. In this section we prove Theo- 
rem 11.8} Corollary 11.91 and Theorem 11.101 

Proof of Theorem 1 1.81 Fix an odd integer k > 3 and let a be a probability measure 
on R with cr({0}) < 1 and J \x\ k da(x) < oo. Denote rrij := J x^da(x) for 1 < 
j < k. If mfc = 0, the theorem is trivial. If < 0, we define a, the "reflection 
through of a", by o(A) := a (—A) for measurable sets A. It is straightforward 
that j x^da(x) = (— l) J mj for 1 < j < k implying that the RHS of the bound 
(J7J) of the theorem is the same for a and a. Since it is also straightforward that 
n%{k) = n®(k), we see that it is sufficient to prove the theorem with a replacing a. 
Noting that j x k da(x) = —rrik > (since k is odd), we shall henceforth assume, 
WLOG, that m k > 0. 

Set a := > (using that m,k-i > since k is odd and cr({0}) < 1). Suppose 

that, for some n, \i := - Y^i=i f° r ( x i)r=i C 1 is a Chebyshev-type quadrature 
formula for a of degree at least k. Letting f(x) := x fc_1 (a — x), we then have that 

(27) / f(x)dp(x) = I f(x)da(x) = am k _i — m k = 0. 
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We continue by noting that, since k is odd, /(0) = and f(x) > for all x G 
(— oo,a) \ {0}. Thus, ( 127]) implies that either = <5 or /i([a, oo)) > 0. However, 
the former option is impossible since a ^ Sq, which implies m 2 > 0, and \x has 
the same second moment as a. It follows that if we denote £ := max{xj}™ =1 , then 
£ > a > 0. Denoting now g(x) := \=r, we have that 

g(x)da(x) = = — 

However, since g(x) > for all x (using that k is odd) and g(£) > g(a) = 1 (using 
that £ > a > 0), it follows that 

- -g(a) < -g(£) < I g(x)d/i(x) = 



n n n J m k 



m k ~ 1 



whence n > — | — as required. □ 



"k-i 



Proof of Corollary \ 1.91 Let da(x) := l[ 0jOO )(x) exp(— x)dx be the exponential dis- 
tribution. Recall that o~([x, oo)) = exp(— x) for x > and that J x^da(x) = j\ for 
j G N. Fix an odd integer k > 3 and define a new measure o~' k by 

<t£(A) := v(yln [0,2fc]) 

for measurable sets A, where c k '■— (1 — exp(— 2/c)) -1 is chosen so that o' k is a 
probability distribution. Define also the rescaling, a k , of a' k to the interval [0, 1] 

by 

a k (A) := a' k (A * 2/c) 

for measurable sets A, where A*2k := {2kx \ x G A}. Noting that a k is absolutely 
continuous, supported on [0, 1] and having density bounded above by Ck for some 
C > 0, we claim that if k is sufficiently large, a k satisfies the corollary (for that 
k). To see this, we shall prove below that if k is sufficiently large then 

From these inequalities we deduce, using that c k — > 1 as k — > oo and that by 
Stirling's formula, (k — 1)! ~ \/2ixk (^— ^) as k — > oo, 

(Jx k da k ) k - 1 _ (Jx k da' k ) k ~ 1 > k k ~ l IT /e\ k 1 (t^ k 



(/ x k ^da k ) k (f x k - x do' k ) k ~ c k 2 k -\k-\)\ V nk \2J ' 2^fk V2 
as k — > oo. Thus, by Theorem 11.81 for sufficiently large k, 



n ° f k )> _J_ f£V 



2%/fc V2, 

as required. It remains only to prove ( 1281 . The second inequality of ( 1281 follows 
from the fact that J" x^da(x) — j\. To see the first inequality, note first that 

(29) / x^da k (x) = c k I j\ — / x^da(x] 

J V J 2k 
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Second, note that for x >2k and j < k we have x J < (2k) 3 exp ( a ~ 2fc ), which can 
be seen by taking logarithms and differentiating. Thus, for j = k — 1 and j = k, 
if k is sufficiently large, 

/■°° , , , . (2k) j f°° / x\ , 2(2A;V j! 

L x rfor(x) - V" L exp (-2) rfx - "V- - 1 

which, when plugged into (|29|) . proves the first inequality. □ 

Proof of Theorem \l.l(A The bounds (fTUl) and (fTTl) follow directly from the bounds 
(jH) and (jUJ) and the definitions of d n (o"o), n° (A;) and n ao (k). The bound (jSJ) follows 
from Bernstein's Theorem [L2] by replicating the Chebyshev-type quadrature given 
by the upper bound of that theorem to each of the two intervals in the support 
of do. To see we will need the following definition and theorems. We call a 
Chebyshev-type quadrature formula (T5]) symmetric if the measures - Y17=i anc ^ 
n ^^=1 are e Q uai - We will use a special case of a Theorem of Forster and 
Ostermeyer [5]. 

Theorem 2.10. ($\ Section 3, Corollaries 1 and 2] If a probability measure a on 
a bounded interval has a density w(x) satisfying w(x) = w(—x) for all x, then for 
each neN, there exists a symmetric Chebyshev-type quadrature having exactly n 
(not necessarily distinct) nodes and degree of accuracy d n (a). 

We will also use a special case of a Theorem of Peherstorfer [20J . 

Theorem 2.11. (]20j special case of Theorem 3.1] There exists C > such that 
for each n G N, if a Chebyshev-type quadrature §2§ for cq has n (not necessarily 
distinct) nodes (xj)" =1 satisfying G (xi)™ =1 , then its degree of accuracy k satisfies 
k<C\n(Cn). 

We remark that the proof of Theorem 12 . 1 II proceeds by taking a polynomial T(x) 
of degree at most k which satisfies |T(x)| < 1 on the support of do and which is 
positive and grows very fast off the support of o~ (a variant of the Chebyshev 
polynomial may be used). Using the facts that the integral of T with respect to 
do is at most 1, that this integral must equal the integral of T with respect to 
the given quadrature formula and that each node of the quadrature formula has 
weight -, one deduces that n must be very large, to offset the contribution of the 
quadrature nodes outside the support of cfq. 

Now fix an odd integer n G N. By Theorem I2.10[ there exists a symmetric 
Chebyshev-type quadrature (j2J) for <To having exactly n nodes and algebraic degree 
of accuracy d n . Since n is odd, the symmetry implies that is one of the nodes 
of this formula. This implies, by Theorem 12.111 that d n (cr Q ) < C\n(Cn) for some 
C > 0, proving ©. □ 

Remark 2.5. As a final remark for this section, we note that it is possible to have 
a sequence of absolutely continuous distributions with n° ak (k) rising as quickly as 
we want with k. However, the densities of these distributions will have very large 
essential supremums. For example, for k = 2m — 1, we can take a distribution 
with m atoms and with the leftmost atom as small as we want. By Bernstein's 
theorem \1.6[ any Chebyshev-type quadrature for it of degree at least k will have at 



SIMPLE UNIVERSAL BOUNDS FOR CHEBYSHEV-TYPE QUADRATURES 



23 



least as many nodes as one over that atom (since the distribution and its Gaussian 
quadrature coincide in this case). Now, we can convolve this distribution with a 
smooth function which is very close to a delta measure to obtain an absolutely 
continuous distribution whose Gaussian quadrature is as close as we want to the 
atomic measure (in the weak topology), so that Bernstein's theorem implies the 
result. 

2.3. Random Chebyshev-type cubatures on the cube. In this section, we 
prove Lemma [1.111 and Theorem 11.121 

Proof of Lemma Suppose that for some n, k > 1 and all e > we have 
Pn,fc,e(c) > 0. If k = 1, a Chebyshev-type cubature always exists for a (placing 
all nodes on the mean of a). Assume k > 2 and fix a sequence Ej — » 0. Since 
Pn,k,e 3 > we can find a measure 

1 - 

•t=i 

such that || Mfc(cr) — M^a^W^ < Sj. These measures must have a converging 
subsequence as j — > oo (in the sense that the location of the atoms converges) to 
some a' := - J^ILi since if any of the atoms goes to infinity we necessarily have 
||-^fc(°j')||oo — °° since k > 2 and each atom carries a fixed weight -. a' is the 
required cubature. □ 

We proceed to prove Theorem 11.121 Recalling the statement of the theorem, 
we first observe that since the Mj are IID vectors in RPoiyDim(M) } t h e central 
limit theorem g i"V6S t licit S n converges weakly to a N(0, E) RV for some matrix 
E. To prove the proposition we would like to show that E is positive definite and 
that a local limit theorem also holds. This will imply that for large enough n, 
the density of S n exists and is uniformly close to that of iV(0,E), whence it is 
uniformly positive in a neighborhood of the origin. 

For a random variable X e M m we write X : M m — > C for the characteristic 
function /(A) := Ke lX ' x . We use the following local limit theorem from [3]. 

Theorem 2.12. Th. 19.1, Ch. 4}) Let (X n ) n >i be a sequence of IID random 
vectors in M. m with KX\ = and positive definite E := Cov(Xi). Let Q n : = 
-^=(Xi + • • ■ X n ), then the following are equivalent: 

(1) Qx G L p (R m ) for some 1 < p < oo. 

(2) For every sufficiently large n, Q n has a density q n and 

lim sup \q n (x) - 4>oM x )\ = 

where 0o,s is the density of a N(0, E) random vector. 

In our case we take Xi := Mj — EMi and we will show that 

(30) M x E L p (R Po ^ Dim ^) for some 1 < p < oo. 

Note that to use the above theorem it may seem necessary to separately show that 
E := Cov(Mi) is positive definite, but this also follows from ( 130]) since if Cov(Mi) 
were singular then X x would be supported in a linear subspace and ([HP]) would 
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not hold, since in that case Xi(p + A) would equal Xi(fJL) for every A orthogonal 
to that linear subspace. 

Hence Theorem II. 121 will follow by verifying ( |30|) . Such estimates are standard 
in the theory of oscillatory integrals but since we could not find this exact result, 
we prove it using standard methods from the book [23] by Stein. Following that 
book, we use the next estimate of Van der Corput to prove what we need. 

Proposition 2.13. ( ]25| Prop. 2, Ch. VIII]^ Suppose </> : (a,5) — > M is smooth 
and satisfies \(f>^\p)\ > 1 for p G (a, b). Then 



< Cj\- 1/j 



when j > 2 or when j = 1 and <$' is monotonic. The bound Cj is independent of 
cf), X, a and b. 

For the case j = 1 we will not be able to ensure monotonicity, so we will use 
instead: 

Lemma 2.14. Suppose (j) : (a, b) — > R is smooth with \4>'(p)\ > 1 for p G (a,b). 
then 



(31) 



b 



iX <t>(p)d 



A A pe(a,fe) 



Proof. The proof is a slight variation on the proof of the previous proposition for 
the case k — 1, as it appears in [25]. Using integration by parts, 

L <w(p) iW(pr L d P \i\4/(p)) 

The boundary terms are majorized by | and the second term satisfies 

For u G let D u denote the directional derivative operator in the direction 
u, and let D{ be its j-th power; i.e., 

d j 

Di(f)(x) = —f(x + pu)\ p=0 . 
We continue with two simple technical lemmas: 

Lemma 2.15. Let Q : R d ->• R be a non-zero polynomial of degree j, then there 
exists u G such that D J U (Q) is a non-zero constant function. 

Proof. Denote m := f'^ 1 ) and let Pf : R d R m be defined by Pf(x) := 

(x a )\ a \=j, where a is a multi-index. We first note that the image of S^" 1 under P d 
is not contained in any proper linear subspace of lR m . This follows since otherwise 
there would exist 77 G S m_1 such that 77 • Pj(u) = for all u G S d_1 contradicting 
the fact that 77 • Pj is a non-zero homogeneous polynomial. 

Now decompose Q as Q = Q\ + Q2 where Qi is a non-zero homogeneous poly- 
nomial of degree j and Q 2 is of degree at most j — 1. Write Qi(x) = J2\ a \=j a a% a - 
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It follows from the above that we may choose u G such that Pf{u) is not 
orthogonal to (a a )\ a \ = j. Hence taking p G R, we see that Q(pu) is a non-zero poly- 
nomial of degree j in p, from whence it follows that for every x G R d , Q(x + pu) is 
a polynomial of degree j in p with the same leading coefficient. Finally, we deduce 
that j^-Q(x + pu)\ p= o is a non-zero constant function as required. □ 

Lemma 2.16. There is c k > such that for every direction rj G §P ol yDim(M)-i ^ e 
function rj ■ P£ satisfies that there exists 1 < j < k and u G § d_1 with 

mm \Dl(r]-P^(x)\>c k . 

xe[-i,i] d 

Proof. Fix r] G SiW^M)- 1 anc i denote P(x) := r? • P^(x). Since P(x) is a 
non-zero polynomial of some degree j < fc, by the previous lemma, there ex- 
ists a u G § d_1 such that -D{-P is a non-zero constant. Hence, in particular, 
mhXj.g^^d |D^(P)(a;)| > 0. The lemma follows since 

max max min {D^rj ■ P^)(x)\ 
ne §d-i i<j<k ie[-i,i]i 

is a continuous function of rj and § Pol yDim(M)-i j s a compact set. □ 

Proof of TheoremEJM Denote /(A) := Mi (A). Fix a direction r? G S Pd y Dim ( fc ' d )- 1 , 
let r > and consider 

(32) f{rri) = Ee ir7rMl = [ e ir7rP * {x) dx. 

J[-i,i] d 

Our goal is to prove an estimate of the form 

(33) \f(rrj)\ <C k r- s 

for some C k and s > independent of rj and r. Such an estimate will imply (|30|) 

for p > PolyDim(M) _ 

Applying Lemma [2. 161 we obtain a number 1 < j < k and a direction w G § d_1 
such that 

(34) min \Di(rj ■ P^)(x)\ > c k , 

where c k > is independent of r/. In addition, we may define 

(35) C k := max max max \D 2 a (rj ■ P k )(x)\ < oo. 

^ g §PolyDim(fc,d)-l j-g^l^jd ue §d-l 

We now decompose our space to the line H = {pu} pe ^ and if -1 . We say that 
y G H L is contributing if there exists p G R such that pu + y G [—1, For 
contributing y's define 

a y = mm{p eR\ pu + y E [-1, 

&„ = max{p G R | pu + y G [-1, 

For non-contributing y's set a y = b y = 0. Note that by a simple I2 estimate, if 
\y\ > y/d, then y is non-contributing. We note that we may estimate the integral 
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321) we are after as 



JrvPi{pu+y) dpdy 



< 



dy. 



1 H- L n[-Vd,Vd] d 

Using the inequalities ( 134")) and (|35|) we may apply Proposition 12.131 and Lemma 
EH to (j)y(p) := (rj ■ Pi){(M + y) to obtain 



(36) 



r<t>y{p) d 



for Ck independent of r], r and y. Plugging this estimate in the previous integral 
we finally obtain 

\f{rrj)\<C k r- l ' k 

as required. □ 

2.4. Local cubature formulas. In this section we prove Theorems 11.131 and 
11.141 The idea behind the proof of Theorem II. 131 is to present the sphere in spher- 
ical coordinates. Partition the spherical coordinate space into suitable boxes and 
then use the fact that the measure on each box is a product measure to construct 
Chebyshev-type quadratures for them using our one-dimensional quadrature re- 
sults. The same idea with a few variations was used in Wagner [27] . To prove 
Theorem 11.141 we use the fact that the measure on the cylinder is a product of 
the measures on the x-axis and the measure on the sphere. We then partition 
the x-axis to small intervals and construct a product Chebyshev-type cubature 
on each interval using our construction for the sphere and our one- dimensional 
quadrature results. 

2.4.1. Sphere Construction. We begin the proof by introducing spherical coordi- 



nates. Let Ang d := (</>, 6 X , . . . , 9 d -i) G R d and tt d := {0 < 
vr for 1 < % < d - 1}. Then define T : Q d ->■ R d+1 (formally T d 

d-i 

sin(0) JJsin(0i), 



P < 2tt,0 < 6i < 

by 



^(Ang,)! 
^(Ang d ) 2 



i=i 

d-l 



COS 



sin 



i=i 



T(Ang d ) 



d-l 

cos(%_ 2 ) Y[ 

i=j-l 



sm 



for 3 < j < d + 1. 



This is a continuous and onto mapping of Q d to S d . Further endowing Q d with 
the measure 



dfj, d (Ang d 



d-i 

n 

i=l 



sm 



l )de l 



the map becomes measure preserving (§ d is endowed with the surface area mea- 
sure o~ d ). We will embed Q d into fl d +i and write (with slight abuse of notation) 



Ang, 



d+l 



(Ang d ,6 d ). Note also that d/j, d+ i( Ang 



d+l, 



sin (9 d )d9 d dfi d (Ang d 
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Similarly, we will embed S d = T{Q d ) into E> d+l = T(Q d+l ) by T(Ang d+1 ) = 
(sin(0 d )T(Ang d ),cos(0 d )). 

We now construct the partition we shall use in Theorem 11.131 In the spherical 
coordinates space f2 d , the sets of the partition will be taken as boxes, that is, 
Cartesian products of intervals. 

Proposition 2.17. For each d > 1 and < r < 1 there exist C(d),c(d) > 
(independent of t), K = K(r,d) > 0, and a partition of Qd (up to measure 0) 
into boxes Dx, . . . , with side lengths smaller than 1, Diam(T(Di)) < C{d)r 
and Hd(Di) > c(d)r d for all i. 

Proof. We proceed by induction. For d = 1, we partition Qi = [0, 2tt] into 
intervals of length 27r/[^ L ] < 1. It is straightforward to see that the required 
properties hold. Assume that the proposition holds for dimension d — 1. We will 
construct boxes Di, . . . , Djc satisfying the required properties for dimension d. 
First partition [0,tt] into m := [-] length — < 1 intervals (/i)™^ 1 (overlapping in 
their end points) by It := [a$, a i+ x) with a« := i^-. For each < i < m — lwe 
will define a set Cj of boxes of the form D x I i where D C Q^-i is a box. Then 
Di, . . . , Dk will be the union of all of the Cj. 

Fix < i < m — 1 and define r := sin( a,+ ^ I+1 ) and r' := min(-, |). We then 
have for any a > 0, 

max sin(6 l ) < Cr, 

ai<8<a i+ i 

(37) [ 1+1 sm. a (e)de > c(a)r a T, 

J at 

ct < rr' < t 

for some C, c > independent of all other parameters and c(a) > depending only 
on a. Indeed, the first inequality follows from the facts that < a, < a i+1 < 7r 
and sin is non-negative, concave and continuously differentiable on this interval. 
The second follows from these facts and |aj + i — | = ^- > > |r. Finally, the 

T 

RHS of the third inequality follows directly from the definition of r' and the LHS 
relies on the facts that sin a; > |x for < x < 1 and \a i+ i — > |r to deduce 
that r > ||r, from which the inequality follows. 

Using the induction hypothesis, let Di, . . . , D/> be the partition of which 
satisfies the proposition for r'. The set C, is the set (Dj x Ii)f =1 . Fix 1 < j < if 
and let £) := Dj x ij. It remains to check that Diam(T(D)) < C(d)r and fJ*d{D) > 
c(d)r d for some C(d),c(d) > independent of r. To check the former, note that 
T(D) = {(sm(9d-i)T(Dj), cos(9 d -i)) \ a, < 6 d -i < cti+i} and hence by the triangle 
inequality, the induction hypothesis and ( |371) . we have 

Diam(T( J D)) < max sin(0 d _i)Diam(T(Z)j)) + cos(aj) — cos(flj_|_i) < 

< Cr ■ C7(d - l)r' + Cr < C(d)r. 
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To check the second bound, note that by the product structure of the measure 
and (13"T|) . we have 



Hd{D) = n±. x (Dj) 
as required. 



sm 



(0 d _i)d0 d _i > c(d-l)(T') d -cid-iy-W > c(d)r 



□ 



For the subsets {Ei} of Theorem 11.131 we take Ei := T(D{) where {Di}f =l are 
the boxes of Proposition 12.171 (with the same d and r as in the theorem). For the 
rest of the proof fix 1 < % < K and, for brevity, denote D := Di and E := E^. Let 
h : M. d+1 — > R be defined by h(z) := z a for a multi-index a with \a\ < k. Since 
D is a box, we may write D := J x I± x • • • x (these {/,} are different from 
those used in the proof of Proposition I2.17p . Note that 



(3? 



h(z)da d (z) 



/i(T(Ang,))^(Ang d ) 



D 



d-1 



' sin Ql (0) cos" 2 (0)d0 TT sin E 

D „_i 



'=^(0 9 )cos 



sin 9 (0,)d0 9 



9=1 



d— 1 



COS 



"9+2 [ 



sin Ql (0)cos O2 (0)#]J / sin 9 

q=l 

We begin the construction of our cubature formula by constructing quadratures 
for the intervals J and I n . 



Lemma 2.18. Given < 7 < 1, integers k>l, 0<q<d — 1 and any 
< <7i < o"2 < T g; where tq := 2tt and r 9 := 7r /or 1 < g < <i — 1, and such that 
a% — a\ < 1, let m > 1 6e t/ie minimal integer such that (^rx) m+ < \- Then there 
exists C = C(d) such that for any integer n > C m there exist (wj)™ =1 Q (01,02) 
satisfying 



7 



(39) 7 f T — — I sn,'' l+ ''(0)<W''-'(0W0 - - > Mil 4 '(//,) cos''-' (//,) < 
J sin 9 (0)a0 

/or a// integers hi, k 2 > suc/i £/iai ki + k 2 < k. 

We first show how to use the lemma, then give its proof. 

Proof of Theorem \1.13[ Denote I := J. Using the lemma, for each < q < d — 1, 
let (wg,jj)™ =1 be the (?/j) satisfying (139]) for the given k and for 7 := Let 



( x i,j)l=i be the Cartesian product (y t 



0,i,j)j=l 



x (y d -i,ij)^ =1 . Finally let 



be defined by z it j := T(xij). Note that for any /i as in 
and using that | sin(0)| < 1 and | cos(6 l )| < 1 that 



•;3>j = l 



we have by ([55]), (1591 



1 

/ h{z)da d {z) - — 
7e n i=i 



< ^7- 



To finish the proof of the theorem, it remains to show that the Zij provide an 
approximate Chebyshev-type quadrature also for g(z) of the form g(z) = (z — w) a 
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for w G S d and a multi-index a with \a\ < k. Fix such a g. For a multi-index 
G (N U {0}) d+1 we write /3 ^ a if /3 g < ct g for all g. Then * 



0"d 



n 



v JJt/ 0^a g=l v ^ 9/ j=l /3^a q=l yt q/ 



< 



q=l 

which completes the proof of the theorem. □ 

Proof of Lemma \2.18l Fix 7 > 0, integers k > 1, < q < d — 1 and fci,A;2 > 
such that k\ + k2 < fc. Fix also < <j\ < 02 < r q satisfying 02 — 01 < 1. Let 
/(0) := sin fcl (#) cos fc2 (fl) and write 



f{9) = <a{0 - ax)' + r m {9){9 ~ °x) m+1 

i=0 

the Taylor expansion with remainder term of f(9) up to degree to. Recall that 

j"(m + l) /M -' 

^m(^) = ( m+1 )i where for > 01, 9 is some number in [ax, 9). By the Cauchy 
estimates, we have for any p > 0, 

M0)| < p~ [m+l) max + 9)\ < p- {m+1) max | e *( fc i+^)(*+0)| < ^-(m+i)^ 

\z\=p \z\=p 
8£(a lt e) 6e(<ri,0) 

Choosing p = 222 ^- we obtain |r m (0)| < (^fi) m+1 - We thus choose to > 1 to be 
minimal such that (^rj-) m+1 < If we now find (%)" = i satisfying 

I 71 1 rO"2 

(40) - J2 Vs = r> • ,, flWfl / sin'Wde 

71 ~t iaa sin* (0)d0 

for all integer < s < to, then these will satisfy the requirements of the 
lemma since 



1 



JJ 2 sin" (9)d0 



f(9)^{9)d9-^Y^f{ yj 



< 



< 



ro* * , fl I" r m {9){9 - ax) m+1 sm q (9)d9 - - ^r m ( % )( % - - ax)"* 1 
J CT1 sm q {9)d6 J ai n ^ 



< 



7, 



where we used \a 2 — o~i\ < 1 and sin q (9) > for G [0,r 9 ] in the last in- 
equality. To find 's satisfying (H0|) . first scale the problem from [01,02] to 
the [0, 1] interval. Theorem 11.41 then shows that (yj)^ =1 exist for any integer 
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n > 75e 4 mM(12eM) (m 1} where M = j^g^p maxgg [aiitJ2] sin 9 (fl). The lemma 
follows since M < C(g) independently of o\ and 02. D 

2.4.2. Cylinder Construction. We prove Theorem II. 141 in the special case W = 1. 
The general case follows from this as follows. If we want the general case with 
parameters k, L', W, r' and 5' we can take the W = 1 construction with the same 
k and parameters L := ^t,W = 1, t := (5 := ^,^ fc and rescale its result by a 
factor of W. 



Proof of Theorem for W = 1 . We first use Theorem 11.131 with its input pa- 



rameters d, k, t, 5 taken to be, in terms of the input parameters to Theorem 11.141 
d — 2, k, r, j^f, respectively. We thus obtain sets Ei, . . . , E K i C E> d ~ 2 and, for 

each 1 < i < K', (zij)j =i C Ei, satisfying the assertions of Theorem 11.131 Next, 
for each 1 < i < X', we define intervals (ii, g ))2i C [—§!/, §L] by the following 
procedure: := a^J with a^o := — with subsequent a^g's defined by 

the rule 

T d-x 



(41) / v(x)dx-- 

J[cH, q -i,a i>q ] 0- d _ 2 {&i) 

and with the integer set to be the maximal one for which [a^ mi -ii a i,m^ Q 
[— |L, |L]. Here we recall that, by Property (1) of Theorem 11.131 we have that 
a ^ d ^ E ^ < C'(d)r for some C'{d) > and, thus, using that 1 < v(x) < 2, r < 1 
and our assumption that L > C(d), we can, and do, take C(d) sufficiently large 
to ensure that m; is well-defined and ai }jn . > L for all i. Note also that the (Ii q ) 
satisfy Diam(/j i9 ) < C(d)r for some C(d) > and all i and q. Finally, we take 

K := Y^f=i m i anc l ^ ne required sets (Dj)f =l to be all sets of the form I^ q x Ei for 
1 < i < K' and 1 < q < mj. 

We continue by establishing properties (I) and (II). By the fact that the (Ei) are 
disjoint up to surface measure 0, and our construction, it follows that the (Dj) are 
disjoint up to z^i-measure 0. Furthermore, as remarked above, we have ai >mi > L 
for all i and, thus, V2L,\ (Pl,i \ yjf =1 Dj)) = 0, proving (I). In addition, as noted 
above, Diam(Jj g ) < C(d)r for all i and q and thus, by the diameter bound on the 
(Ei) and our construction, it follows that Diam(Z)j) < C(d)r for some C(d) > 
and all j. Relation (T4T]) ensures that V2L,i(Dj) = r d ~ x for all j and this implies 
the bound on K by a volume estimate. Thus, (II) is proved. 

It remains to establish (III). Fix 1 < i < K', 1 < q < rrii and, for brevity, 
denote D := I is x Ei. We start by constructing the points (w 0,3)^=1- We apply 
Theorem 11.41 to the interval Ii A with the weight v restricted to that interval (by 
rescaling the interval to [0, 1] and renormalizing the measure to be a probability 
measure) and obtain (xj i9J )"^ 1 C [a^-i, a^g] satisfying 
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for any integer < r < k. Furthermore, no may be any integer such that no > 
75e 4 JfcM(12eM) ( * -1) where M := X'f^'. 1 max xe[ai _ lfii ^v(x) < 2. Hence, 



•K'o-i < x ) dx 



q ] v [ x ) 



there exists C > 0, independent of all other parameters, such that no may be any 
integer satisfying n > C k . Finally, we take the points (wd j)?=i to be the Carte- 



sian product of (xi t q j j)^l 1 and (zij)j =l (of Theorem II .131) . This implies 



where N 



n 



d-2 



and ni can be any integer satisfying m > C\(d) 



n := noiV, 

m (d-2,k,S/(4L) k ) 



for some d(d) > 0. Note that m (d - 2, k, 5/ (4L) fe ) > k (using that L > C(d) 
and taking C(d) > 1), implying that n can be any integer of the form n^" 1 for 
n 2 > C 2 (rf)™o(d-2,fc,5/(4L) fe ) ; for some > 0) as required. 

It remains only to verify the approximate quadrature condition of (III). Fix 
h : M. d — > R, := (w — y) a for ?/ G P2L,i an d a multi-index a with |a| < fc. 

Let h : M^ 1 — >■ R be defined by h(w) := (w — where y G -R^ -1 is defined by 

:= (so that y G S d_2 ) and a G (N U {O})^ 1 is defined by ct^ := a i+ i. Using 
the product structure of the (wdj), D, v-il,\ and h, and (H2j) and (|14l) . we obtain 



1 " 



i 



h(w)dv 2 L,i(w) 



1 ^ - 



■ 8 

(twi - yi) ai w(wi)ciw;i 



1 



< max 

w,y£P 2 L,i 



wi ~ V\ 



Od-2(Ei 



h(w)da d - 2 (w] 



< 



5 



(4L) k 



<6, 



where we used that L > C(d) and took C(d) > 1 in the last inequality. 



□ 



3. Open questions 

(1) What is the best possible upper bound on n° a {k) and n a {k) using only the 
information contained in R a or using other simple properties of a? Can 
the conclusion of Theorem 11.31 be improved? In our lower bounds section, 
we showed that Theorem 11.31 is sharp, up to constants, for measures with 
bounded density. Is this also the case for other measures? 

(2) For which measures a does n®(k) or n a {k) grow only polynomially with 
kl Can such behavior be deduced using only simple properties of a (i.e., 
without knowing the orthogonal polynomials of <r)? Is it always the case 
for measures on a finite interval, with a density which is bounded away 
from and infinity? 

(3) Salkauskas [23] puts a probability measure on the set of length k mo- 
ment vectors (normalized volume measure) from which he deduces that 
the proportion of those vectors for which a Chebyshev quadrature (i.e., a 
Chebyshev-type quadrature with n nodes and algebraic degree at least n) 
exists is exponentially small in n. Can Salkauskas' result be extended to 
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give the typical degree of a n node Chebyshev-type quadrature? Is this 
typical degree a power of n or logarithmic in n? 

(4) Obtain quantitative theorems for random Chebyshev-type quadratures. 
Do they help to show existence of Chebyshev-type quadratures? In par- 
ticular, for the case of the uniform measure on the cube, what is the order 
of magnitude of Nq in Theorem 11.121 ? (see also remarks preceding the 
theorem) 

(5) Theorem 11.31 and its extensions give upper bounds for n a {k) for measures 
a supported on a finite interval. Can similar theorems be written for 
measures in higher dimensions? or for measures with unbounded support? 

(6) Theorems 11.31 and 11.51 give upper bounds for n a {k) in the cases when a 
has at least k + 4 atoms, or a non- atomic part. However, we know from 
Theorem II. II that Chebyshev-type quadratures exist once we have roughly 
| atoms. Can we bound n a (k) using simple properties of a for measures 
having between | and k + 4 atoms? 
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